Finger-gate array quantum pumps: pumping characteristics and mechanisms 
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We study the pumping effects, in both the adiabatic and nonadiabatic regimes, of a pair of finite 
finger-gate array (FGA) on a narrow channel. Connection between the pumping characteristics and 
associated mechanisms is established. The pumping potential is generated by ac biasing the FGA 
pair. For a single pair (N — 1) of finger gates (FG's), the pumping mechanism is due to the coherent 
inelastic scattering of the traversing electron to its subband threshold. For a pair of FGA with pair 
number N > 2, the dominant pumping mechanism becomes that of the time-dependent Bragg 
reflection. The contribution of the time-dependent Bragg reflection to the pumping is enabled by 
breaking the symmetry in the electron transmission when the pumping potential is of a predominant 
propagating type. This propagating wave condition can be achieved both by an appropriate choice 
of the FGA pair configuration and by the monitoring of a phase difference (j> between the ac biases 
in the FGA pair. The robustness of such a pumping mechanism is demonstrated by considering a 
FGA pair with only pair number N = 4. 



I. INTRODUCTION 

Quantum charge pumping (QCP) 

has become an active field in recent 
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concerned with the generation of net transport of 
charges across an unbiased mesoscopic structure by 
cyclic deformation of two structure parameters. Original 
proposal of QCP, in the adiabatic regime, was due 
to Thoulessi^ He considered the current generated 
by a slowly varying travelling wave in an isolated 
one-dimensional system. The number of electrons 
transported per period was found to be quantized if 
the Fermi energy lies in a gap of the spectrum of the 
instantaneous Hamiltonian. Aiming at this quantized 
pumped charge nature of the adiabatic pumping, Niu 
proposed various one-dimensional periodic potentials for 
the adiabatic quantum pumping (AQP)^ and pointed 
out the importance of the quantized charge pumping in 
utilizing it for a direct-current standard^ 

Another way to achieve the AQP was suggested by 
Hckking and Nazarovj 3 - who studied the role of inelas- 
tic scattering in the quantum pumping of a double- 
oscillating barrier in a one-dimensional system. Intended 
to stay in the adiabatic regime, they invoked a semi- 
classical approximation and had assumed that the Fermi 
energy ep 3> hQ, where is the pumping frequency. 
This semiclassical treatment of the inelastic scattering is 
known to be inappropriate for the regime when either the 
initial or the final states are in the vicinity of the energy 
band edge. Such a regime, however, is our major focus 
in this work. It is because the coherent inelastic scatter- 
ing becomes resonant when the traversing electron can 
make transitions to its subband threshold by emitting 
mhflml£l Depend on the system configuration, this and 
another resonant inelastic scatterings will be shown to 
dominate the pumping characteristics^ 

A recent experimental confirmation of AQP has been 



reported by Switkes et alS Two metal gates that defined 
the shape of an open quantum dot were ac biased^ 3 , with 
voltages of the same frequency but differed by a tun- 
able phase differenced DC response across the source and 
drain electrodes is the signature of the AQP. This has 
prompted further intensive studies on AQP in: quantum 
dotSfSiS double-barrier quantum wells^ pumped volt- 
age?ii noiseless AQPji^ heat current, 13 incoherent pro- 
cessesr^i^ quantum rings^^^ and interacting wiresi^ 

An alternate experimental effort in generating AQP 
involves surface acoustic wave (SAW)i 24 i 25 i 26 i 27 i 28 Gener- 
ated by an interdigitated SAW transducer located deep 
on an end-region of a narrow channel, the SAW prop- 
agates to the other end-region of the narrow channel 
while inducing a wave of electrostatic potential inside the 
channel. Electrons trapped in the potential minima are 
thus transported along the narrow channel. Both Mott- 
Hubbard electron-electron repulsion in each such trap 
and the adiabaticity in the transport are needed to give 
rise to quantization in the pumped current^ As such, 
the channel has to be operated in the pinch-off regime^ 

In this work, we propose to study yet another experi- 
mental configuration for QCP in a narrow channel. The 
proposed configuration consists of a pair of finite FGA, 
with the number N of FG's in each FGA being kept to a 
small number. In contrast to the SAW configuration, the 
FGA pair sits on top of the narrow channel, rather than 
locating at a distance far away from it; and the most sig- 
nificant QCP occurs in regimes other than the pinch-off 
regime. The FG's orient transversely and line up longi- 
tudinally with respect to the narrow constriction. As is 
shown in Fig. ^ pumping potential can be generated by 
ac biasing the FGA pairs with the same frequency but 
maintaining a phase difference 4> between them. Since 
the wave of electrostatic potential induced in the narrow 
channel is directly from the FG's, rather than via the 
SAW, our proposed structure has the obvious advantage 
that the working frequency is not restricted to the fre- 
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FIG. 1: (color online). Top view of the proposed system 
structure for the case of pair number N — 4. A FGA pair 
locates on top of a narrow channel. Vi denotes the amplitude 
of the potential energy, and <j> is the phase difference. 



quency of the SAW, usg = 2nvs/d. Here vg is the phase 
velocity of the SAW, and d is the pitch in the FGA. Fur- 
thermore, when the working frequency is different from 
tos, the contribution from SAW to the pumped current 
will be negligible. 

Below we shall show how the ac biased FGA pair plays 
a subtle role in the generation of QCP. In Sec. II, we 
present our theoretical model for the FGA pair calcula- 
tion of the pumped current generated by the FGA pair 
configuration. In Sec. Ill, we present the pumping char- 
acteristics and demonstrate that resonant coherent in- 
elastic scatterings are the underlying pumping mecha- 
nisms. Finally, in Sec. IV, we present our discussion and 
summary. 



II. FGA PAIR MODEL 

The potential V(x,t) in a narrow constriction induced 
by a FGA pair is represented by 

JV 

V(x, t)=Yl V ^ cos ( m ) + V ^ x ) cos (^ + 0)> W 

where N is the number of FG's per FGA. We assume 
that the ac biased FGA pair are localized, respectively, 
at positions Xi and Xi+Sxi, namely that Vu(x) = V\S(x— 
Xi) and V2i(x) = V2<5(x — x% — Sx) with a relative phase 
difference 4>. These FG's are evenly spaced, with a pitch 
d, and are located at xi — (i—l)d for one FGA and Xi+Sx 
for the other. The relative shift between the FGA pair 
is Sx = a d, where the fractional shift < a < 1 . In the 
following, we consider the case of the same modulation 
amplitude V\ = V2 = Vq. Depending on the choice of the 
values for <fi and a, V(x, t) will either be predominantly of 
a propagating or a standing wave type. A sensible choice 
can be made from considering the lowest order Fourier 



component of V(x,t), given by 

Vi = —r- {cos AT 2; cos ilJ + cos [K{x — Sx)] cos(f2J + (/))}, 

" ^ ( 2 ) 

where K = 2n/d. For our purposes in this work, an 

optimal choice is (j) = ir/2 and a = 1/4, in which V(x, t) 

is a predominant left-going wave. 

The Hamiltonian of the system is H = H y + H x (t), in 
which H y = —d 2 /dy 2 + u^y 2 contains a transverse con- 
finement, leading to subband energies e n — (2n + l)u) y . 
The time-dependent part of the Hamiltonian H x (t) is of 
the dimensionless form H x (t) — —d 2 /dx 2 + V(x,t). Here 
appropriate units have been used such that all physical 
quantities presented are in dimensionless form^ 

In the QCP regime, the chemical potential /i is the 
same in all reservoirs. Thus the pumped current, at zero 
temperature, can be expressed asifi 

I = ~ [ dE[T^(E)-T^(E)}. (3) 

Here the total current transmission coefficients include 
the contributions by electrons with incident energy E in 
incident subband n, which may absorb or emit mVL to 
energy E m = E + mQ by the FG pumping potentials, 
given by 

Afe— 1 00 

T -(-)( £; )=E E T n ^ } (E m ,E), (4) 
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where Afs stands for the number of occupied subbands. 
The summations are over all the propagating components 
of the transmitted electrons, and includes both the sub- 
band index n and the sideband index m. The subscripted 
arrow in the total current transmission coefficient indi- 
cates the incident direction. These coefficients are calcu- 
lated numerically by a time-dependent scattering-matrix 
method-iS*^ 



III. NUMERICAL RESULTS 

In this section we present the numerical results for 
the pumping characteristics of cither a single FG pair 
(JV = 1) or a finite FGA pair (JV > 2). In these two 
cases that the pumping characteristics are due to dif- 
ferent resonant inelastic scattering processes. For def- 
initeness, the unit scales in our numerical results are 
taken from the GaAs-Al x Gai_ x As based heterostructure. 
The values that we choose for our configuration param- 
eters are u> v — 0.007, subband level spacing Ae = 2oj y 
(~ 0.13 meV), d = 40 (~ 0.32 /im), and V = 0.04 (~ 
28.7meV A). From the value of Vq, and the assumed FG 
width ~ 0.05 /im, the amplitude of the potential induced 
by a FG is - 0.057 mV. 
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A. single FG pair case 

In this subsection we investigate the pumping charac- 
teristics for the case of a single FG pair. FigureElpresents 
the dependence of the total current transmission coeffi- 
cients on the incident electron energy \i. We replace the 
chemical potential \x by 



1 



X » ~ As ' 2 ' 



(5) 



which integral value corresponds to the number of propa- 
gating subbands TV's in the narrow channel. The pumping 
frequency is higher in Fig.^a), with £1 = 0.6Ae (f2/27r ~ 
18 GHz), than that in Fig. Efb), where ft = 0.1 Ae 
(U/2tt ~ 3 GHz). We select the phase shift <f) = ir/2 
and a = 1/4. 
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beneath a subband bottomiSi The peak structures in 
of Fig. Hfb), and at A M = A/"s + 0.2, are 2f2 resonant 
structures. 




FIG. 3: The pumped currents versus with the same pa- 
rameters used in Fig. The solid, and dashed, curves corre- 
spond, respectively, to Q — 0.6Ae and Q, — 0.1 Ae. 



In Fig. |3 we can see that T^(X fl ) does not equal to 
T—,(X I1 ), this allows the occurrence of the pumped cur- 
rent. Moreover, between integral X^ values, T<_ > T_> 
on the left region of a dip structure, while T<__ < on 
the right region of the dip structure. This has an impor- 
tant bearing on the dependence of the pumped current 
on fi, as is shown in Fig. [3J The pumped current rises, 
and drops, on the left, and right, region of a X^ip, re- 
spectively, in accordance with the relative changes in T_» 
and about the same X,a p . Hence the peaks of the 
pumped current depend on the pumping frequency, at 
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FIG. 2: Total current transmission coefficient versus X M for a 
pair of FG; (a) Q = 0.6Ae and (b) Q = 0.1 Ae. The transmis- 
sion of the right-going (left-going) electrons are represented 
by the solid (dotted) curve. The subband level spacing is Ae. 
Parameters a — 1/4 and <j> — n/2 are chosen to meet the 
optimal condition. 

At integral values of X^, the total current transmission 
coefficients (^(X^) exhibit abrupt changes. This is 
due to the changes in the number of propagating sub- 
bands in the narrow channel. Between integral X^ val- 
ues, T—> («_) both show dip structures. The dip struc- 
tures are located at Adi P = A/"s + 0.6 in Fig. [21 a), and at 
-^dip = A/s + 0.1 in Figs.EJb). These dip structure loca- 
tions are the same for both T_> and T<_, and are resonant 
structures associated with inelastic scattering that causes 
an electron to jump into a quasibound state (QBS) just 



reassuring us that the pumping is dominated by the 
aforementioned resonant inelastic process. 

Besides the trend that the pumped current in Fig. 
drops with the pumping frequency, we would like to 
remark on a more interesting result: that both the 
adiabatic and nonadiabatic behaviors can be found in 
the same curve. Since the adiabatic condition is given 
by S> f2, the curve for Q, = 0.1 Ae in the regions 
M s + fl/Ae < X^ < J\f s + 1 corresponds to the adia- 
batic regimes, while the other X^ regions are nonadia- 
batic regimes. This is checked also with our other calcu- 
lation, which is not shown here, using the Brouwer ex- 
pression^ For the higher pumping frequency, f2 = 0.6Ae, 
the adiabatic condition is not satisfied in the entire X^ 
region, even though the pumping characteristics resem- 
ble that of the adiabatic one in the regions Af s + il/Ae < 

Xfj, < M S + 1. 
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B. finite FGA case 

In this subsection we present the numerical results for 
the pumping characteristics of a finite FGA pair. QCP 
for two prominent modes of tuning the system are con- 
sidered. These are (i) tuning of the electron density by 
the back-gate technique, and (ii) tuning of the channel 
width by split-gate technique. 

1. tuning back-gate 

We present the numerical results for the pumping char- 
acteristics of a FGA pair with N = 4 that is realized by 
the back-gate technique. The dependence of the total cur- 
rent transmission coefficients on A M is shown in Fig.^J in 
which the pumping frequencies are (a) SI = 0.6Ae, and 
(b) Q = O.lAe. The choice of the parameters d, (f>, and a 
is the same as in the previous subsection, but the latter 
two parameters give rise here to an equivalent left-going 
wave in the pumping potential V(x,t). 
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FIG. 4: Total current transmission coefficient versus X M for 
N = 4; (a) SI = 0.6Ae and (b) SI = 0.1 Ae. The transmission 
of the right-going (left-going) electrons are represented by the 
solid (dotted) curve. The parameters a = 1/4 and cj> = n/2. 

The curves in Fig. 0] show additional structures, other 
than the dip structures that has been discussed in the last 
subsection. These additional structures are valley struc- 
tures that occur at different values for T_^(A M ) and 



T<_. In a region between two integral values of A M , the 
valley structure of T^{X I1 ) occurs at a lower A M . This 
shows clearly the breaking of the transmission symme- 
try by the pumping potential. Furthermore, the valleys 
are separated by AX M = £1/Ae. This can be understood 
from resonant coupling conditions Ek = £k-K — and 
£k+K = £k — fi for, respectively, the right-going and the 
left-going k. From these conditions, the valley locations 
are at 
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where the upper sign is for positive, or right-going, k. 
These locations, expressed in terms of X^, are given by 

and are at X^ = 1.19, 1.79, 2.19, 2.79, 3.19, and 3.79 
for the case of Fig. Ha), and A M = 1.39, 1.49, 2.39, 2.49, 
3.39, and 3.49 for the case of Fig.QjT). The matching be- 
tween these numbers and our numerical results in Fig. 0] 
is remarkable. In addition, energy gaps open up at these 
k± locations, causing the drop in the transmission and 
the formation of the valley structures^ All these results 
reassure us that the time-dependent Bragg's reflection is 
the dominant resonant inelastic scattering in our FGA 
pair structure. 

On the other hand, the adiabatic condition is here 
given by e gap 3> £1, where £ gap is the effective energy gap 
of the instantaneous Hamiltonian^ Since e gap is given by 
the widths of the valley structures, therefore contribu- 
tions of the valleys to the pumped current is nonadiabatic 
in Figs. Ola), because the valleys are well separated, and 
adiabatic in Figs.0Jb), because the valleys overlap. 
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FIG. 5: Pumped current versus X M . The choices of parame- 
ters are the same as in Fig. 0] The solid, and dashed, curves 
correspond, respectively, to SI = 0.6Ae and SI = O.lAe. 

In Fig. [5] we present the A M dependence of the pumped 
current for the cases in Fig.^J The pumped current peaks 
at X^ that lies in the middle between a valley in T_^(X M ) 
and the corresponding valley in TU(A^). The locations 
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which depend on both the pitch d and the pumping fre- 
quency fl. The peaks have flat tops for the solid curve, 
when Q = 0.6Ae. Comparing with the total current 
transmission curves in Fig. Ufa), we see that the flat- 
topped peak profile is associated with the complete sep- 
aration between the valleys in T_> and T«_ ■ This is in the 
nonadiabatic regime. In contrast, for the case when the 
valleys overlap, such as in Fig.0Jb), the pumped current 
peaks no longer carry a flat-top profile, as is shown by the 
dashed curve in Fig. [SJ This is in the adiabatic regime. 
Meanwhile, their peak values are lowered, ft is because 
cancellation sets in when the valleys overlap. We note 
that the pumped currents are of order nA. 

The robustness of the time-dependent Bragg reflection, 
on the other hand, is demonstrated most convincingly 
by the number of charge pumped per cycle at the max- 
imum /Max of the pumped current. In the dashed curve 
of Fig. [SI the pumped charge per cycle per spin state 
Q P = (27r/0)/ M ax/2e = 0.495, where J M ax = 0.48 nA 
and f2 = 0.1 Ae = 3.03 GHz. To get a unity, or quan- 
tized, charge pumped per cycle per spin state, one can 
fix the pumping frequency £1 = O.lAe, N = 4, (j> = n/2, 
and d = 40, then tune the other pumping parameters 
as V Q = 0.09 and a = 0.15 to obtain Q P = 0.992 at 
= 3.465 (not shown here). In this frequency regime, 
the pumping would be expected to be adiabatic, accord- 
ing to Thoulessi and Niu£ when e gap 3> J7. However, in 
our case here, the energy gap is at best only partially 
opened, as we can see from the nonzero transmission in 
Fig. Efb), because we have only N = 4 FG pairs. Thus 
our result shows that the condition of occurrence of the 
of AP is less stringent than we would have expected orig- 
inally. 2 In other words, the pumping effect of our FAG 
configuration is robust. 

It is also worth pointing out that the pumped currents 
are positive in Fig. showing that the net number flux 
of the pumped electrons is from right to left. This is con- 
sistent with the propagation direction of the electrostatic 
wave in V(x,t)Mi 



2. tuning split-gate 

Thus far, we have explored the dependence of the FAG 
pair's QCP characteristics on by the use of the back- 
gate technique. Another way of tuning the QCP charac- 
teristics is via the modulation of the channel width (or 
subband level spacing Ae). This can be realized experi- 
mentally by the use of the so-called split-gate technique. 
Hence we present, in Fig. [S] the transverse confinement 
dependence of both the total current transmission coeffi- 
cients and the pumped current. The transverse confinc- 
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FIG. 6: The dependence on subband level spacing Ae of (a) 
the total current transmission coefficient, and (b) the pumped 
current. The abscissa is depicted by Eq. I|1O0 where fi = 
0.049 and N = 4. Pumping frequency Q = 0.0084 in all 
curves except for the dotted curve in (b), where fl = 0.0014. 
Parameters cf> = tt/2 and a = 1/4 for all curves except for 
the dashed curve in (b), where a = 1/5. In (a), the solid 
(dashed) curve is for T_»(X g ) [T_(X S )], and contributions 
from the second Fourier component of V(x,t) are indicated 
by arrows. 



ment is depicted by 
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which is linearly related to the effective channel width, 
and that its integral value corresponds to the number of 
propagating subbands in the channel. In this mode of 
tuning the QCP characteristics, [i is kept fixed. 

In Fig. [ffia). except for /i, which is fixed at 0.049, 
and Lo y , which varies with X g , other parameters such 
as fl = 0.0084, 4> — 7r/2, and a — 1/4 are the same as in 
Fig. Ha). The solid (dashed) curve is for (T<_). Both 
the BS and the time-dependent Bragg reflection features 
are found. The expected locations of the BS, given by 
the expression 



X„ 



li 



1 



2J /x-O' 



(11) 



are at 1.1, 2.3, and 3.5, and they match the BS loca- 
tions in Fig.EJa) perfectly. Here n is the subband index. 
The expected locations of the valleys, associated with the 



6 



time-dependent Bragg reflection, are given by the expres- 
sion 

^=H n+ £)^4 ! (i2) 

thus they should be at X g = 1.03, 2.1, 3.14 for T^(X g ), 
and at X g = 1.15, 2.4, 3.73 for T^(X g ). Again, they 
match the valley locations in Fig. Efa) remarkably. 

Besides, there are in Fig. EJa) two additional valley 
structures, indicated by arrows, at which T_>(X g ) and 
T^(Xg) fall one on top of the other. These structures do 
not contribute to the pumped current, and they are due 
to the time-dependent Bragg reflection from the second 
order Fourier component of V(x,t). The second Fourier 
component of V(x,t) is in the form of a standing wave, 
given by cos (2 Kx) [cos fit + sin fit]. That both of the 
additional valleys all appear in T^(X g ) and T^(X g ) can 
be understood from the fact that more resonant coupling 
conditions come into play for the case of standing wave. 
The resonant coupling conditions are Ek — su±2K ± ^, 
and £fc = £k±2K T ^- As such, the valley locations are 
given by the expression 

x ^ = l\ 1+ v Ji r\ (13) 

for n = 0, and for e± = [K(l T fl/(2K) 2 )} 2 . Accord- 
ingly, these 2K time-dependent Bragg reflection valley 
locations are expected to be at 1.36 and 1.73, which co- 
incide with the two additional valleys in Fig.^a), which 
are indicated by arrows. We note, in passing, that con- 
tributions from higher Fourier components diminishes, as 
is seen by comparing the valleys from the first and the 
second Fourier components of V(x,t). 

The X g dependence of the pumped current for the case 
in Fig.HJJa) is represented by the solid curve in Fig.HJb). 
The peaks have flat tops because the valleys in the cor- 
responding T^(X g ), T^(Xg) are well separated. The 
pumped current for £1 = 0.0014, the same frequency 
as in the case of Fig. 0Jb), is depicted by the dotted 
curve in Fig. EJb). The peaks are not flat-topped and 
the magnitudes are much smaller because the transmis- 
sion valleys overlap. For comparison, we also present the 
case when parameter values differ slightly from that of 
the optimal choice. As is shown by the dashed curve in 
Fig. EJb) , where all parameters are the same as for the 
solid curve except that a is changed from 1/4 to 1/5, 
the basic pumped current peaks in the solid curve re- 
main intact. This demonstrates the robustness of the 
QCP against the deviation in values of the configuration 
parameters from the optimal choice. 

Interestingly, there are two additional features in the 
dashed curve of Fig. EJb) : namely, an additional pumped 
current peak at X g — 1.5, and an increase in the peak 
value for the pumped current near X g = 3.5. That 
both of these features are found to arise from the second 
Fourier component of V(x,t) is supported by the out- 
come of our analysis performed upon the Fourier compo- 
nent of V(x, t). This method of analysis has thus far been 



successful in providing us insights on the pumping char- 
acteristics presented in this work. The m-TtH Fourier 
component of V(x,t), apart from a constant factor, is 
given by the form 

V m = {[cos(TO7ra) — sin(m7ra)] cos[mKx' — fit — 7r/4] 
+ [cos(m7ra) + sin(m7ra)] cos[mKx' + fit + 7r/4]} , 

(14) 

where x' = x — 6x/2. V m consists, in general, of waves 
propagating in both left and right directions. But when 
a = 1/4, as we have discussed before, Vi becomes a pure 
left-going wave and V2 becomes a pure standing wave. 
The case of a = 1/5, however, have both Vi and V2 
consisting of waves in opposite propagation directions. 
Therefore, in contrast with the a = 1/4 result, additional 
contributions from the 2K Bragg reflection are expected 
for the case a = 1/5. This additional contribution should 
peak at the mid-point between two transmission valleys 
for the 2K Bragg reflections, and the expression for X g 
is given by 

*»=1+Hfe (15) 

where tu = K 2 + (fl/2K) 2 . For the case of the dashed 
curve in Fig. HJb), the values of X g = 1.54 and 3.6 are 
shown to match the locations of the additional features 
nicely. Finally, we can extract information of the sensitiv- 
ity of the pumped current characteristics to a by looking 
at the coefficients of the left-going and right-going waves 
in V m . For a = 1/5, the coefficients of V\ for, respec- 
tively, the right-going and the left-going waves are 0.22 
and 1.4. This shows that V\ is still dominated by the 
left-going wave and thus explains the tiny modifications 
to the pumped current peaks at X g = 1.1, and 2.3. But 
for V2, the coefficients for, respectively, the right-going 
and the left-going waves are —0.95 and 1.57. This shows 
that V2 deviates quite significantly from that of a stand- 
ing wave, and so explains that the additional peaks from 
the 2K Bragg reflections are quite large. 

IV. DISCUSSION AND SUMMARY 

It is interesting to note in passing that our proposal of 
the FGA pair configuration is different, in three aspects, 
from the voltage lead pattern proposed earlier by Niu£ 
First of all, the pumping mechanisms to which the con- 
figurations are catering to are different. It is the mech- 
anism of translating the Wannier functions in a given 
Bloch band in Ref. Q, while it is the mechanism of the 
time-dependent Bragg reflection in this work. The former 
mechanism is adiabatic by nature but the latter mecha- 
nism is shown, in this work, to hold in both the adiabatic 
and non- adiabatic regimes. 

Second, the configurations are different in the number 
of sets of voltage leads invoked. A third set of voltage 
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leads was instituted by Niu to fix the Fermi energy at the 
middle of the instantaneous energy gap in order to main- 
tain the adiabaticity of the pumping. Since our interest 
here is on the general pumping characteristics, including, 
in particular, their dependence on the Fermi energy, it 
suffices us to consider a simpler configuration — the FGA 
pair configuration. Third, the number of voltage lead 
expected, and needed, in a voltage lead set is different. 
Our results demonstrate the resonant nature of the time- 
dependent Bragg reflection, and that the pumping char- 
acteristic is robust — requiring only a FGA pair with small 
N. Hence the FGA pair configuration proposed in this 
work should be more accessible experimentally. 

In conclusion, we have proposed a finger-gate array 



pair configuration for the generation of quantum charge 
pumping. Detail pumping characteristics have been ana- 
lyzed, the robustness of the time-dependent Bragg reflec- 
tion in QCP has been demonstrated, and the pumping 
mechanism is understood. 
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